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Abstract. In this paper, we introduce a new notion named as Schrodinger soliton. So- 
called Schrodinger solitons are defined as a class of special solutions to the Schrodinger 
flow equation from a Riemannian manifold or a Lorentzian manifold M into a Kahler 
manifold N. If the target manifold N admits a Killing potential, then the Schrodinger 
soliton is just a harmonic map with potential from M into N. Especially, if the domain 
manifold is a Lorentzian manifold, the Schrodinger soliton is a wave map with potential 
into N. Then we apply the geometric energy method to this wave map system, and 
obtain the local well-poscdness of the corresponding Cauchy problem as well as global 
existence in 1 + 1 dimension. As an application, we obtain the existence of Schrodinger 
soliton of the hyperbolic Ishimori system. 



1. Introduction 

In this paper we intend to study a class of special solutions of the Schrodinger flows 
from a Riemannian manifold or a Lorentzian manifold into a Kahler manifold. First, let 
us recall some preliminaries on Schrodinger flows. Let (M, g) be a Riemannian manifold 
or a Lorenzian manifold and (N, h, J) be a Kahler manifold, where J denotes the complex 
structure and h is the Kahler metric. The Schrodinger flow is a map w : R x M — y N 
which satisfies the equation 

(i.i) ( If = JHrW ' 

[ w{0) = w . 

where r(w) = trace g 'V 2 w is the tension field of w, and wq is an initial map from M to N. 

The Schrodinger flow from a Riemannian manifold stems from fluid mechanics and 
physics. It is a problem with strong physical backgrounds and a long history. A century 
ago Italian mathematician Da Rios studied the motion behavior of vortex filament and 
discovered the well-known Da Rios equation which can be formulated as 

It = Is x 7 SS , 

where 7(s,t) : S 1 x R — y R 3 is a closed space curve for a fixed time t. By differentiating 
the above equation with respect to s we obtain the so called ferromagnetic spin chain 
system which is just the Schrodinger flow into S 2 . For the existence theory of Schrodinger 
flow from a Riemannian manifold, we refer to [1, 5, 6, 7, 10, 28, 29, 33] and references 
therein. Yet, for the Schrodinger flow from Lorentzian manifolds, little is known. In 1984, 
Ishimori [18] proposed a model as a 2 dimensional analogue of the classic continuous 
isotropic Heisenberg spin chain, which also describes the evolution of a system of static 
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spin vortices in the plane. The hyperbolic-elliptic Ishimori problem is a spin field model 
with the form: 



(1.2) 



d t s = s x Ds + b(d x s ■ d y 4> + d y s ■ d x 4>), 
A(ft = 2s ■ (d x s x d y s), 

where s : M 2 x R ->• S 2 ^ E 3 , □ = d 2 x - d 2 , b G R and lim^^i^ s(x, y, t) = (0, 0, -1). 
The Cauchy problem associated to Ishimori system (1.2) has been studied extensively 
in the past decades, see for example [22, 30] and references therein. When 6 = 0, this 
system gives a simple example of Schrodinger flow from a Lorentzian manifold. 

Kenig, Ponce and Vega [23] have ever studied the following Schrodinger equation which 
is analogous to the Schrodinger flow from Lorentzian manifold: 

{Ou 
— = i££u + P(u, Vu, u, Vu), 
u(0) = u . 

where u = u(t,x) is a complex valued function from R x W l , .if is a non-degenerate 
second-order operator 

i<k j>k 

for some k G {1, • • • , n}, and P : C 2n+2 — y C is a polynomial satisfying certain constraints. 
They proved the local well-posedness of the above initial value problem in appropriate 
Sobolev spaces. 

Since it is difficult to establish a general existence theory for Schrodinger flow from 
Lorentzian manifolds, we return to looking for some special solutions. We recall that 
in [11] the authors proposed to study the periodic solutions of the Schrodinger flow in 
the case where the target manifold iV is a Kahler-Einstein manifold with positive scalar 
curvature. If the target manifold is just the standard sphere S 2 , they employed the well- 
known symmetric variational principle to show the existence of some special periodic 
solutions to the flow from a closed base surface with convolution symmetry. In particular, 
they needed to reduce the Schrodinger flow to a elliptic equation and established the 
following lemma on reduction. 

Reduction Lemma. Assume there exists a non-trivial holomorphic Killing vector field 

V on N, and let S t be the one-parameter group of holomorphic isometries generated by 

V with So = I, the identity map. Then w(t) = S t o u with u : M — > iV is a solution to 
(1.1) if and only if u is a solution to the equation 

(1.4) t{u) = -J{u)V{u). 

Proof. Directly computing by the definition of tension field, we get 

t(w) = r(St ou) — dS t o t(u) + r(S t )(du, du). 

Since S t is an isomorphism, we have r(S t ) = and hence 

t(w) = dS t o t{u). 

On the other hand, 

d 

w t = Qj:(St ou) = V(S t ou) = dS t o V(u). 

The last equality holds because the single parameter group St satisfies S t o S s — S t + S - 
Differentiating this at s = 0, we get dS t o V — V(S t ). 
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Next, because V is holomorphic, i.e. [J, VV] = 0, we have 

J o dS t = dS t o J. 
Combining above equalities together, we arrive at 

(1.5) w t = dS t o V(u) = J(w)t(w) = J (St o u)dS t o r(u) = dS t o J(u)t(u). 

dSt is an isomorphism on the tangent space, so (1.5) is equivalent to (1.4). □ 

It is easy to see that the special solution to Schrodinger flow given by the above lemma 
is some kind of solitary wave solution. In fact, for a linear Schrodinger equation defined 
on a flat torus T m 

iw t = Aw, 

a solitary wave solution is of the form w = ue lkt where k is a positive constant, and v is 
a real function which satisfies the equation Au + ku = 0. Here, e lkt can be viewed as a 
holomorphic isometric group with one parameter. Therefore, we define the Schrodinger 
soliton as follows 

Definition. Suppose u is a solution to (1.4) derived in the Reduction Lemma, then w(t) = 
S t o u is called a Schrodinger soliton solution of (1.1). 

A solution to equation (1.4) is a map with prescribed tension field. In general it is hard 
to solve the equation because the elliptic system is not of a variational structure. There 
are only a few results under some strong assumptions, see [4] for example. 

However, if there exists a smooth function A e C°°(N) on N , such that JV = VA is 
the gradient vector field of A, then the equation becomes 

(1.6) t(u) = -VA(u), 

and it's easy to see that this equation is the Euler-Lagrange equation of the following 
functional: 

(1.7) F(u) = E(u) - I A(u)dV g . 

J M 

Here 

E(u) = \f \Vu\ 2 dV g 
is the energy functional of maps u e W 1 ' 2 (M,N), where \Vu\ = trace g (u*h). 

In the case M is a Riemannian manifold, the solutions to equation (1.6) are harmonic 
maps with potential A from M into N. Once we have the above variational structure, 
many powerful tools which are adopted to study harmonic maps work for the present 
problem and many results on harmonic maps can be extended. For formal results on 
harmonic maps with potential, we refer to [2, 3, 12, 13]. 

In this paper, however, we focus on the situation where the base manifold is Lorentzian. 
It is well-known that the hyperbolic harmonic maps from a Lorentzian manifold are usu- 
ally called wave maps and the well-posedness of wave maps has been intensively studied 
by many mathematicians; see for example [31, 34, 35] and many references therein. We 
will see below that the Schrodinger soliton from a Lorentzian manifold (or Lorentzian 
Schrodinger soliton for short) satisfies a perturbed wave map equation. It's worthy 
pointing out that this kind of wave map with potential emerges naturally as a simplified 
equation of the dynamics of weak ferromagnets magnetization when N = S 2 [16]. 

Indeed, let (M\,gi) be a compact Riemannian manifold with the Riemannian metric 
9\ = gafidx a dxP and M = Rx be a Lorentzian manifold equipped with a Lorentzian 
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metric g = dt 2 —g±. Denote the covariant derivative for functions on Mi and M by V and 
V respectively. We will always embed the compact target manifold N into a Euclidean 
space ~R K . Then the equation (1.4) becomes 

(1.8) Du = A(u)(Vu, Vu) - J(u)V(u), 

where □ = d 2 — A is the wave operator, Vu — u t + Vu and A(u)(-, •) is the second 
fundamental form of iV C WL K . Using the Christoffel symbols T k j of N, one can write 
explicitly in local coordinates that 

(A(u)(Vu, Vu)) k = Y%u\u\ - g^V^Vpui. 

Equation (1.8) is a nonlinear wave system. In particular, if there exists a Killing potential 
(See Section 2 for the definition) A G C°°(N) such that JV = VA, the equation becomes 

(1.9) Du = A(u)(Vu,Vu) - VA(u). 

We will call a solution to equation (1.9) a wave map with potential. We will consider 
initial data 

(1.10) (u(0),ut(0)) = (u ,«i); u x (x) G T uo(x) N, a.e. x G M 1 

and study the corresponding Cauchy problem. Our main result is the following theorem: 

Theorem 1.1. Suppose (Mi,gi) is an m- dimensional compact Riemannian manifold 
with m > 1 and M = Ix Mi is equipped with a Lorentzian metric g = dt 2 — gi, let 
N be a compact Kdhler manifold with a Killing potential A such that VA = JV. Let 
k > m = [y] + 1, where [y] denotes the integer part of y. Than for initial maps 
u G W k ' 2 (Mi,N) andui G W k ~ l > 2 (Mi, T U0 N), the Cauchy problem (1.9), (1.10) has 
a unique local solution u satisfying u G L°°([0, T), W k ' 2 (Mi, N)) and u t G L°°([0,T), 
W k ~ 1,2 (Mi, T U N)). Moreover, if the initial data is smooth, so is the solution. 

Remark 1.1. Although for the sake of consistency with the Schrddinger soliton, we only 
discuss wave maps with Killing potentials in this paper, by exactly the same procedure 
one can verify that Theorem 1.1 holds for wave maps with any potential A, i.e. for any 
smooth function A : TV — >■ R. 

In the classical wave map theory, it has been shown that the Cauchy problem of 
wave map is locally well-posed on Minkowski space Rxl" with initial data {uq,Ui) G 
W k ' 2 (W l , N) x W k - 1 ' 2 (R n ,T U0 N), where k = ^ for n > 3 and k > § for n = 2. (See 
Theorem 7.2 in [31].) On the other hand, the C°°-regularity of wave equations is well- 
known by the theory of paradifferential operators. Thus Theorem 1.1 is a generalization 
of the well-known results for wave maps to the current perturbed wave map system on 
Lorentzian manifolds. Note that mo is the critical exponent on the manifold M, since 
there are no fractional Sobolev spaces on manifolds. 

This generalization won't take much effort since the perturbing term is of lower order. 
However, in this paper, we employ a new method, namely, the geometric energy method 
which first appeared in Ding and Wang's work [10] to tackle this problem. It's worthy 
to point out that the geometric energy method is a powerful tool in dealing with various 
kinds of geometric evolution equations. It's also the first time shown in this paper that the 
wave map (with potential) can be handled by this method. It provides a simplified and 
uniform method which avoids the complicated analysis of fixing moving frames, choosing 
Columb gauge, etc. (See [32] for example.) 

Another advantage of this method is that we can directly obtain the C°°-regularity of 
the solution to the Cauchy problem with smooth initial data. The fact that the Cauchy 
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problem is locally well-posed with initial data in W k ' 2 for all k > m dose not directly 
imply the local well-posedness in C°°. Because the space of smooth maps C°°(M 1 , N) = 
n^ =mo W k ^ 2 (M 1 , N) itself is not a Banach space, and the standard techniques such as fixed 
point theory do not apply here. Our method provides an uniform lower bound of the 
maximal time Tk for all k > mo + 1, see Lemma 3.6 below. With this bound, we are able 
to assert the existence of a local solution u G C°°([0, T) x Mi, N) to the Cauchy problem 
with smooth initial data. Moreover, the maximal time T only depends on the geometry 
of N, II 

M o||iy m o+ 1 ' 2 an d ||wi|| vK m o. 2 , see Theorem 3.7. 
In addition, we prove the global existence of solution to the Cauchy problem (1.9), 

(1.10) on 1 + 1 dimensional Lorentzian manifolds. This is an analogous result to the wave 
map theory, see [14] and [31]. 

Theorem 1.2. Let M\ = S 1 and N be a compact Riemannian manifold. Suppose A is 
a smooth function on N and u G W 2,2 (S l , N) , Ui G W 1 ' 2 ^ 1 , T U0 N), then the Cauchy 
problem (1.9), (1.10) has a unique global solution u G L°°(R + , W 2,2 (S 1 , N)). 

Therefore, in this case by Reduction Lemma, we get a Schrodinger soliton solution 
w(t) — S t o u to the Schrodinger flow (1.1) with initial data w(0) = u, which is a special 
global solution. 

Corollary 1.3. Let M = M.x S 1 be a Lorentzian manifold and N a compact Kahler man- 
ifold with a Killing potential A. Then given u G W 2 ' 2 (S 1 , N) and u\ G W l,2 (S l , T UQ N), 
there exists a Schrodinger soliton solution of (1.1). 

Remark 1.2. Theorem 1.1 also holds for Lorentzian space M = IR n+1 . Similarly, The- 
orem 1.2 and Corollary 1.3 hold true if we replace S 1 by R 1 . One may check this by 
following the argument in [10] . The crucial fact is that the interpolation inequality in 
Theorem 3.1 is scaling invariant, and hence the main estimate in Lemma 3.6 does not 
depend on the diameter of the domain. 

Particularly, since there is an natural Killing potential on S' 2 (see Section 2), Corol- 
lary 1.3 provides a global solution to the Ishimori system (1.2) when 6 = 0. Indeed, if 
S 2 G M 3 is the standard sphere and A(-u) = u 3 denotes the projection of u = (ui, w 2 , u 3) £ 
S 2 to the z-axis in R 3 , then A, which is the first eigenfunction of S 2 , is a Killing potential. 
The gradient field of A is given by 

VA(u) = P(u)e 3 , 

where e 3 = (0, 0, 1) is the unit vector and P(u) denotes the orthogonal projection of M 3 
to T U S 2 . Precisely, 

P(w)e 3 = e 3 - (it, e 3 )u. 
The complex structure of S 2 is J{u) = ux. Thus the corresponding Killing field V is 

V(u) = -J(h)VA(m) = —u x P(u)e 3 = -u x e 3 , 

and the isometry family St is just a rotation around the z-axis. Then if s(t,x,y) = 
S t o u(x,y), by Reduction Lemma, the equation of s in Ishimori system (1.2) can be 
reduced to 

(1.11) Ou = (\d x u\ 2 - \d y u\ 2 )u - P(u)e 3 . 

Thus given initial data (u ,Ui), we get a solution u on Lorentzian space R 1+1 to equa- 
tion (1.11) by Theorem 1.2. Hence s(x, y, t) = S t ou(x, y) gives a global periodic solution 
to the initial value problem of Ishimori system (1.2) with s(0) = u. We call such a solu- 
tion s a Schrodinger soliton solution of Ishimori system (1.2). Therefore, we obtain the 
following corollary. 
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Corollary 1.4. 7/6 = 0, toen awen w e W 2 ' 2 ^ 1 , JV) and Ul e W^flR 1 , T U0 AT) ; i/jere 
exisfc a Schrddinger soliton solution of Ishimori system (1.2). 

The rest of the this paper is organized as follows: in Section 2 we briefly introduce the 
Killing potential; in Section 3 we prove Theorem 1.1; finally we prove Theorem 1.2 and 
hence Corollary 1.3 and 1.4 in Section 4. 

2. Killing potential and some remarks 

We know that the Schrddinger soliton equations are not of variational structure gen- 
erally. So, it is very difficult to solve (1.4), since the classical variational methods can 
not be used to approach this problem. In fact, it may do not admit any solution at all. 
Then a natural question is: when dose the equation (1-4) have a variational structure? 
One has found the question relates closely to whether a Kahler manifold admits a Killing 
potential function or not. Therefore, let's recall the notion of Killing potential as follows. 

Definition. If A is a smooth function on a Kahler manifold (N,J), and the gradient 
field of A has the form: 

VA = JV, 

where V is a Killing field on N, then A is called a Killing potential. 

Obviously, if there exists a Killing potential on (N,J), then (1.4) is of the desired 
variational structure. Now, a question confronting us is what kind of manifolds do admit 
Killing potentials? Fortunately, one has made great progress on the existence of Killing 
potentials on a Kahler manifold in differential geometric field. Recently, Derdzinski and 
Maschler studied the so-called special Kahler-Ricci potentials which is a special kind of 
Killing potential, and gave a local classification for the Kahler manifolds admitting such 
potentials. It's also related to the conformally-Eintein Kahler metrics. One can refer to 
[8, 9, 19] for more details. 

For completeness, here we give several basic lemmas about Killing potential. 

Lemma 2.1. ([9]) Suppose A is a smooth function on a Kahler manifold, then the follow- 
ing conditions are equivalent: i) A is a Killing potential; ii) VA is a holomorphic vector 
field; Hi) V 2 A is Hermitian. 

Proof. Let V = — JVA, then A is a Killing potential is equivalent to say V is a Killing 
potential, which means VV is skew-symmetric, i.e. 

(2.1) (W)* + W = 0. 

Since W = -JV 2 A, (V 2 A)* = V 2 A and J* = -J, (2.1) is equivalent to 

V 2 A oj -Jo V 2 A = [V 2 A, J] = 0, 

which means VA is holomorphic. Thus i) and ii) are equivalent. On the other hand, if 
V 2 A is Hermitian, i.e. V 2 A(X, JY) = -V 2 A( JX, Y) for any vector fields X, Y. Then 

V 2 A(X, JY) = (X,V y JVA) = -V 2 A(JX,F) = -(V X JVA,V Y >. 

This is equivalent to the skew-symmetry of V = — JVA, which is equivalent to i). □ 

Lemma 2.2. ([9]) Suppose (N,h,J) is a Kahler manifold. If Hi(N,M) = 0, then for 
every holomorphic Killing field V there exists a Killing potential A, such that VA = JV . 
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Proof. Since V is Killing and holomorphic, VV is skew symmetric and commutes with J. 
Thus if we let W = JV, then VW is symmetric. This implies the corresponding 1-form 
£ = L W h is closed, since 

(dO(X, Y) = h(V x W, Y) - h(X, V Y W) 

for any vector fields X, Y. So there exist a function A such that dA = £ and hence 
VA = W = JV. □ 

In fact, the existence of Killing potential is a complicated problem and somehow re- 
lated to the topology of the underlying manifold. The following lemma gives a sufficient 
condition for the existence of Killing potential: 

Lemma 2.3. Let A be a C°° funcion on a Kdhler manifold (M,g) such that 
(2.2) V 2 A + X Ric = o-g, 

where Ric is the Ricci tensor, and x, & ar & some C°° functions. Then A is a Killing 
potential. 

Proof. It is a direct corollary from iii) of lemma 2.1 and the fact that Ric and g are 
Hermitian. □ 

From this lemma, one can see that there are plenty of manifolds admitting Killing 
potentials, including special cases of independent interest. For example, compact Kahler 
manifolds with function A satisfying (2.2) for constants Xi a sucn that x° > are known 
as Kahler- Ricci solitons ([27], [39]). Also, Riemannian manifolds admitting functions A 
satisfying (2.2) with x — have been studied extensively, and their local structure is 
completely understood in [21]. 

We know that it is always an important issue that how many closed geodesies exist 
on a compact Riemannian manifold. An one-dimensional Schrodinger solitons from S 1 
into a compact Kahler manifold with a Killing potential A is a geodesic with potential. 
Since A is closely relevant to the geometry and topology of the target manifold, it is of 
significance that we study the existence of such geodesies. Naturally, we may ask the 
following 

Question 1: At least how many closed geodesies with potential A exist on 
a closed Kahler manifold with Killing potential? 

On the other hand, we should mention another important special case. When N is 
a compact Kahler-Einstein manifold with positive scalar curvature, it is known that for 
every Killing field V, JV = VAi is the gradient vector field of the first eigenfunction 
Ai of the Laplace-Beltrami operator Ajy on iV ([20]). By virtue of this fact and Sacks- 
Uhlenbeck's perturbed technique, Ding and Yin [11] proved there exists an infinite number 
of inequivalent periodic solutions to the Schrodinger flow (periodic Schrodinger solitons) 
from S 2 into S 2 (see also [17]). In this case the potential function in the above Question 1 
is just the first eigenfunction on N. In fact, more generally we may consider the following 

Question 2: Let N be a closed Riemannian manifold and Ai(x) be the first 
eigenfunction of the Laplace-Beltrami operator A N . At least how many 
closed geodesies with potential Ai(x) exist on N? 

3. Local well-posedness 

In this section, we will use the geometric energy method in [10] to prove the local well- 
posedness of Lorentzian Schrodinger solitons into a compact Kahler manifolds with a 
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Killing potential and wave maps with potential. We need to recall an important theorem 
proved in [10]. This is a generalized Gagliardo-Nirenberg inequality. 

Let 7r : E — > Mi be a Riemannian vector bundle over an m-dimensional Riemannian 
manifold Mi and let D denote the covariant derivative on E induced by the Riemannian 
metric. Then we can define a Sobolev norm via the bundle metric for every section 
s e T(E) by 



k 



\\ s \\H k 'i = il-^ s IU«- 

1=0 

Theorem 3.1. ([10]) Suppose s e C°°(E) is a section where E is a vector bundle on M±. 
Then we have 

(3- 1 ) II-^IIlx - C ll S llff*.« (MIl'"' 

where 1 < p, q, r < oo, and j/k < a < l(j / k < a < 1 if q — m/(k — j) ^ 1) are numbers 
such that 

1 j 1 A 1 k, 

- = - + - + « ■ 

p m r q r m 

The constant C only depends on Mi and the numbers j, k, q, r, a. 

For Lorenzian manifold M = Mx M x with metric g = dt 2 —gi and the compact manifold 
iV which is embedded into ~R K , let D denote the covariant derivative on the pull-back 
tangent bundle u*(TN) over Mi of u e C°°(Mi, N) and D = D t + D denote the covariant 
derivative on the bundle over M. Recall we also use V and V to denote the covariant 
derivative of functions on M x and M respectively. For convenience we denote Du = Vw 
and Du = Vw. Obviously, D 2 u = (V 2 w) T is the tangent part of V 2 w. 

Then by the theorem, for Du € T(u*(TN)), we have 

(3.2) II^^ILp < C W Du \\h^ P M llI" a • 

Ding and Wang also showed that the H h ' p norm of section Du is equivalent to the normal 
Sobolev W k+l,p norm of the map u. Precisely, we have 

Lemma 3.2. ([10]) Assume that k > m/2. Then there exists a constant C = C(N,k) 
such that for all u G C°°(Mi, N), 

k 

ii^^ii w k ~ x * 2 — ^ ^ ii-^^iigfc- 1 . 2 
i=i 

and 



k 

" \W 



\\Du\\ Hk - lfl < CJ2 II Vl 
i=i 

Now we return to the equation (1.4), using the covariant derivative D, we can rewrite 
the equation: 

m 

(3.3) t(u) = trace fl ( J D 2 u) = D 2 t u - D a D a u = -J(u)V{u). 

a=l 

To prove the existence of the above equation, usually one needs to choose a suitable 
approximate equation for which the existence is easy to prove, and some uniform a priori 
estimates of solutions with respect to the parameter e needs to be established. Here we 
follow [40] due to Y. Zhou and use the viscous approximation 

(3.4) D 2 t u - D a D a u - eD a D a u t = -J(u)V(u), 



where e > is a small parameter. Or equivalently, 

(3.5) uu - eAu t -Au + J(u)V{u) = A{u){Vu, Vu) - eT{u){Au t )±T u N, 
where T{u) denotes the orthogonal projection to the normal bundle at u, i.e. 

T(u)(Au t ) — Au t — (Au t ) T . 

We already know that 

(Au t ) T = trace gi D 2 u t 

= tiace gi D(Vu t — A{u){u t )) 

= Au t - A(u)(Vu t , Vu) - div(A(u)(u t , Vu)). 

Thus we have 

(3.6) T(u)(Au t ) = A{u){Vu t , Vu) + dw{A{u){u t , Vu)). 

This equation (3.5) may be viewed as a parabolic system for u t . Indeed, the local exis- 
tence and uniqueness of smooth solutions to (3.5) for initial data (-u ? M i) £ C°°(M 1 ,TN) 
such that 

(3.7) (it,itt)(-,0) = (uo,«i) 

can be derived by a fixed point argument using the heat kernel of Mi (see the appendix). 
Actually, Miiller and Struwe [25] used this approximation method to prove the global 
existence of weak solutions to the wave map equation in 1 + 2 dimensions with finite 
energy data. 

We can define the energy density for a map u : M —¥ N and V£ G K. by 

e(t) :=\\Vu(t)\\ 

where 

\Vu(t)\ 2 = \u t (t)\ 2 + \Vu(t)\ 2 . 

Notice that the norm here is different from the norm induced by the Lorentzian metric 
g = dt 2 — gi. This is a convention in wave map theory which we will adopt through out 
this paper. 

Now we define the energy functional for all maps u G W 1,2 {M, N) and Wt G M. by 

E(t) := / e{t)dV gi) 

J (t}xMi 

For this energy functional, we have the following energy inequality: 

Lemma 3.3. For any e G (0, 1], suppose u G C°°(M 1 x [0,T e ), N) is a local solution to 
Cauchy problem (3.4), (3.7). Then we have 

E(t)<E(0)- [ [ (u t ,J(u)V(u)). 

JO JMx 

Particularly, if JV = VA is the gradient field of a Killing potential A, we have 

(3.8) E(t) < E(0) - [ A(u(t)) + / Hu(0)). 

J Mi J Mi 
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Proof. Using the equation (3.5), we have 
dE(t) 



dt 



= / (u tt ,u t ) + (Vu, Vu*) 

J Mi 

= / (u tt ,u t ) - (Au,u t ) 

= / (eAut - J(u)V(u) + A(«)(V«, Vu) - eT(u)(Au t ),u t ) 

JM X 



= -ef |V« t | 2 - f 

JM 1 JMj. 

<- f (J(u)V(u),u t ). 

J Mi 



(J(u)V(u),u t ) 



Integrating this equality from to t, we get the lemma. 



□ 



Thus given a smooth initial data, we can get a local solution u t G C°°(T t x M 1? N) for 
every e > which satisfies the energy inequality. Next, in order to establish the local 
existence of the equation (3.3), we need to derive some uniform a priori estimates for 
solutions u e with respect to e. For this, we denote for a fixed time t G [0, T e ) 

Du 2 = / (Du,bu)= / (D t u, D t u) + (Du, Du) . 
i 2 (^i) J Ml J Ml 

Note again this norm is not the one induce by the Lorentzian metric. 

In the following we will assume Mi is flat, i.e. the Riemannian curvature of Mi vanishes 
identically, to simplify the computations. For the general case, the additional terms 
involving the curvatures of Mi actually do not provide additional difficulties, since the 
derivatives of u appearing in these terms are of lower orders and the curvature of M x are 
bounded. 

Let a be a multi-index with length |a| 
direction, we compute 

Id ~ 2 
2~dt 



D a Du 



I, and D a be the mult i- derivative of space 
{D a Du, D t D a Du). 



Mi 



(3.9) 

2dt L 2 (Mj) 

Changing order of the covariant differentiation, we have 

(3.10) D t D a Du = D a D t Du + D b R(u)(D c u, D d D t u)D e Du, 

where R is the curvature tensor of iV and the summation is taken for all multi-indexes 
b, c, d, e with possible zero lengths, except that |c| > always holds, such that 

(b, c, d, e) = cr(a) 

is a permutation of a. If we denote the curvature terms like the second term on the right 
hand side of (3.10) by Q, i.e. 



then we have 
1 d 
Ydt 

(3.11) 



D a Du 



2 

£2(Afi) 



Q(X,Y) = J2 D bR{u){D c u,D A X)D e Y, 
(D a D t Du + Qi, D a Du) 



Mi 



[D a D 2 t u,D a D t u) + (D a D t Du, D a Du) + (Q 1 , D a Du), 



Mi 
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where Q\ = Q(D t u, Du). 

For the second term in (3.11), we have 



(D a D t Du, D a Du) 



Mi 



(3.12) 



= / (DD a D t u + Q 2 , D a Du) 

= - [ (D a D t u, DD a Du) + (Q 2 , D a Du) 

J Mi 

= - [ (D a D t u, D a DDu + Q 3 ) + (Q 2 , D a Du) 

J Mi 

= / (D a D t u, D a DDu) - (D a D t u, Q 3 ) + (Q 2 , D a Du) 

J Mi 



where Q 2 = Q(Du, D t u), Q% = Q(Du, Du). 

To simplify the notations, we will put all the curvature terms Qi together and use Q 
to denote the sum of those terms. 

Combining (3.11) and (3.12) together and using the equation (3.4), we get 



1 d 
2di 



D a Du 



< / (D a D 2 t u- D a DDu,D a D t u) + \Q\\D a Du\ 

i 2 (^0 J Mi 



Mi 



{eD a DDu t - D a J(u)V{u),D a D t u) + \Q\\D a Du\ 

(eDDD a u t + eQ 4 + tDQ 5 - J{u)D a V{u), D a D t u) + \Q\\D a Du\ 



Mi 



/ 

J Mi 



—e(DD a D t u, DD a D t u) - (J(u)D a V(u), D a D t u) + \Q\\D a Du\ 



<C I \D a u\\D a D t u\ + \Q\\D a Du\, 

J Mi 

where Q 4 = Q(Du, Dut),Q§ = Q(Du, D t u). Obviously, we have 

IQI < igii + ig 2 | + ig 3 | + e|Q4i + e| J Dg 5 | 

(3.13) < C\Q(Du, Du)\ + e\Q(Du, Du t )\ + e\DQ(Du, u t )\ 
< Cj2\ Djl Du\---\D jb Du\, 

where the summation is over all indexes ■ ■ ■ ,jb) satisfying 

(3.14) 3l > j 2 > ••• > j b , l>ji>0, ji + ---+j b + b<l + 3, b>3. 
Thus, we get 



1 d 
2~dt 

Hence 
(3.15) 
1 d 

2~dt 



D a Du 



L 2 (Mi) 



<C I \D a u\\D a D t u\ + C V I \D l Du\\D jl Du\---\D j 

JMi J Mi 



Du\. 



D l Du 



<C I \D l u\\D l D t u\ +(7V / \D l Du\\D jl Du\---\D jb Du 

J Mi J Mi 



L 2 {Mi) j Ml 
I + 11. 
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For convenience, we denote s = Du. Then we can apply Theorem 3.1 on s which is a 
section of the bundle u{t)*TN on M 1 to get 

(3.16) II^MIlp - C IMIff** IMIiT 1 > 

where 1 < p, q, r < oo and j/k < a < 1 satisfy 

,« , ^ 1 7 1 ,1 1 k . 
3.17 _ = A + _ + a . 

p m r q r m 
Let's first estimate the first term / in (3.15). By Holder inequality, 
(3.18) / < C\\D l u\\ L 2\\D l D t u\\ L 2 < CUD'-^HiallD'sllLa. 

Then using the interpolation inequality (3.16), we have 

II r~il— 1 II - ^ — y 1 1 \\a ii ||1— a 

\\D S\\ L 2 < C \\S\\ H l,2 \\s\\ L 2 , 

where a = (/ — I) /I by (3.17). So we get 
(3-19) I <C\\s\\^' 1 \\s\\^\\D l s\\ L 2. 

Next we treat the second term in (3.15), i.e. 



11 = / \D l s\\D 3l s\---\D ]b s\, 

J Mi 

where the indices satisfy (3.14). Here we directly apply Ding- Wang's lemma in [10]. Let 
m = [y] + 1, where [y] is the integer part of y. 

Lemma 3.4. ^[10]^) If 1 < I < m , there exists a constant C = C(Mi,l) such that 

II <C\\s\\A mo , 2 \\ s \\f 2 \\D l s\\ L 2, 

where A = [/ + 3 + (m/2 — 1)6 — m/2]/m and B = b — A. 

Lemma 3.5. ([10]) If I > m , there exists a constant C = C(Mi,l) such that 
(i)ifji = I, 

^<ci| S |i^: o 2 || S |iir /mo n^iii 2 , 

(a) if h < h 

II < C(l+ \\s\\ 2 H u2){l + ||s||^-i, 2 ), 

where A = A(m, I). 

Now we can prove our main lemma. Note that previous computations do not depend 
on the variational structure. But to get the bound on energy, we need to assume that 
JV = VA in the following context. 

Lemma 3.6. Suppose (u ,Ui) G C°°(M 1 ,TN), then there exists 

T = T(\\Vu \\ Hmo ,2, 11^11^0,2) > 

independent of e G (0, 1], such that if u e G C°°(M 1 x [0, T e ], N) is a solution to (3.4), (3.7), 
then T t > T, and 

(3.20) \\Du\\ Hk ,2 < C(||Vuo||h*.3, ||ui||jyk,2),Vt G [0,T], 

for all k > m Q . 

12 



Proof. We still denote s = Du, then the energy functional in Lemma 3.6 is E(t) = i||s||£2. 
Since A is a smooth function on a compact manifold N, it's bounded. From the energy 
inequality (3.8), we have 

(3.21) ||s|| i2 = 2E(t) < 2E(0) - 2 / A(u(t)) + 2 I A(u(0)) < C. 

Now we turn to (3.15). We first consider the case 1 < I < m . According to (3.19) and 
Lemma 3.4, we have 

l±\\ D i s \\ 2 <I + II 

2 alt 11 WL^Ah)- 

< C MIS 1 \\D l s\\ L2 + \\4h^M\D 1 s\\ L 2. 
Summing this inequality from / = 1 to / = m , we get 

Id ii n2 ^ rvV^ II II ('-!)/' i II n A ( 6 >0 mi ii 



i.e. 



(3-22) — ||s||^m ,2 < C(^2 ll s ll//^o' ) ,2 + ^2 W s \\h^)- 



b,l 



where 
If we let 

we have 



(3.23) < £ C! ^°' 



A(b, I) = [I + 3 + (m/2 - 1)6 - m/2]/m . 
/(*) = ||«||fl™o. a + !> 

d 
it 

[ /(0) = ||VM ||H m 0.2 + ||til|| H m ,2 + 1. 



where 



A Q = max{(l -1) /I, A(b, I)}, 

b,l 



and the constant C only depends on ||Vwo||h"»o. 2 , H^iH^o. 2 and the manifolds M ± ,N. 
It follows from ordinary differential equation theory that there exists 

T = T(||VM ||ff"H>,2, ||Wl|| H m ,2) > 

and a constant K such that f(t) < K, i.e. 

(3.24) \\Du(t)\\ H m 0: 2 <K,Vte [0,T]. 

Next we treat the case k > m . (3.15), (3.18) together with Lemma 3.5 leads to 

j t \\ Dl 4h < c\\s\\ 2 Hk , 2 + cJ2( l + \\4 2 h^)(i + INI^-i,*)- 

Summing up from I — 1, • • • , k, we get 

(3.25) J t ^ Hk ' 2 ~ C H {1 + H s ll^^)(l + ||s||£ fc -i, 2 ). 
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Then we perform a induction for k > m . Specifically, we first consider k = m + f. 
From (3.24), (3.25), we get 

j t \\s\\ 2 Hmo+1 ,i < CKJ2(1 + ||s||^ 0+ i, 2 ),Vt e [0,T]. 

By Gronwall's inequality, we get 

\\Du(t)\\ HmQ+ i,2 < c",vt e [0,T]. 

Then by induction, for any k = m > 1 it follows from (3.24), (3.25) that 

\\Du(t)\\ H K2 < C k ,Vt e [0,T\. 

where C& only depends on || Vu \\ H k,2, H M i||,f/ fc >2. 

Thus we proved the lemma. □ 

Now we can prove the local existence of the solution to Cauchy problem (1.9), (1.10) 
with smooth initial data. 

Theorem 3.7. Suppose (u ,ui) G C°°(M 1 ,TN), then there exists 

T = T(\\Vu \\ H m ,2, \\ Ul \\ Hmo ,2) > 

such that the Cauchy problem (1.9), (1.10) has a local solution u G C°°(Mi x [0,T], N). 

Proof. For any e > 0, there is a smooth solution u e G C 0O (M 1 x [0,T e ],N) to (3.4). 
Moreover, u e satisfies the estimate (3.20) in Lemma 3.6 and there is a constant T > 
such that T t > T, Ve > 0. It follows form Lemma 3.2 that 

(3.26) max ||Z>uJ| W fc,2 < C(||Vuo||H fc > 2 > ||wi||.ff*,2), V/c > m , 
te[o,T] 

where the constant C is independent of e. Thus, by letting e — > and applying Sobolev 
embedding theorems, we can find a limit map u G C°°(M 1 x [0,T],N), such that u e — > 
u in C k (Mi x [0,T],N) for any k. It's easy to verify that u is a smooth solution to 
equation (1.8). □ 

From (3.26), one can easily see that the limit map u also satisfies the same estimate, 
i.e. 

max ||Z>u||wfc,2 < C(||Vwo||,f/fc,2, \\ui II ^,2), V/c > m , 
te[o,T] 

In fact, we can say more about u. Namely, the above inequality also holds for k = m — 1. 

Lemma 3.8. Suppose u is a solution to Cauchy problem (1.9), (1.10) given by Theo- 
rem 3.7, then 

(3.27) max ||.Du|| H fc,2 < C(||VwoI|h-*.2, ll^ill H k A, Vfc > — 1- 
te[o,T] 

Proof. The proof goes almost the same with the proof of Lemma 3.6, except for a more 
refined estimate on the curvature term. The observation is that without the approximat- 
ing term eDDu t , there are only three terms left in the curvature term (3.13). Indeed, 
this term becomes 

|Q| < |Qi| + IQ2I + |Q 3 | < C\Q(Du,Du)\ 
< Cj2\D jl Du\---\D jb Du\, 
where the summation is now over all indexes (ji, • • • , jj,) satisfying 

(3.28) Ji>j2>--->Jb, l>ji>0, ji + ---+j b + b<l + 2, b>3. 
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The key is that the sum of the index in (3.28) is I + 2, which is one order lower than I + 3 
in (3.14). With this change, one can verify that all the estimates in the rest part of proof 
of Lemma 3.6 holds for m — 1 instead of m . □ 

Now we are ready to prove the main theorem. 

Theorem 3.9. Suppose u e W k ' 2 {M 1) N),ui € W k ~ 1 ' 2 (M 1 , T U0 N), where k > m . Then 
the Cauchy problem (1.9), (1.10) has a local solution u E L°°([0, T],W k ' 2 (M 1 , N)) with 
iH e L oo {[0,T},W k - 1 ' 2 {M 1 ,TN)). 

Proof. Since uq € W k ' 2 (Mi, N) with k > mo larger than the borderline m/2 for Sobolev 
imbedding into C°(Mi, N), we can approximate uq by smooth maps in C°°(Mi, N)(see 
[10] for a proof). Namely, we may select a sequence (u % 0) u\) £ C ao (M 1 ,TN), such that 

Uq ->■ u in W k ' 2 (M u N),u\ ->■ mi in W k ~ 1 ' 2 (M 1 ,M 2K ). 

Then for any i > 1 and initial data (mq,i4), there exits a local solution -u* which satisfies 
(3.27). Since as i — > oo 

||MollvK fe - 2 — ll M o||w fe > 2 

ll M lll W k - 1 ' 2 ~~ ^ ll M l|| W k ~ 1 ' 2 i 

the estimate (3.27) is uniform with respect to i and only depends on ||i/ ||w*,2 an d 
||mi|| W k-i,2- Hence 

(3.29) max 11^,2 < C(||-uo||iy=. 2 ) H^iIIm^- 1 . 2 ), 

te[o,T] 

(3.30) max ||uj|| W fc-i,2 < C(||-uo||w/ fc . 2 ) IKillw*- 1 . 2 )- 

Therefore we can find a subsequence which we still denote by u\ such that 

tf^um ^([O^lW^iM^N)), 

u\ -± u t in L oo ([0,T},W k - 1 ' 2 (M 1 ,TN)) 

where — 1 denotes the weak * convergence. 

The limit u is a strong solution to (1.8). To show this we only have to verify that for 
any v G C°°{M 1 x [0,T],M x ), there holds 

(3.31) f f (Du- A(u)(Vu,Vu),v) = - I [ (J(u)V(u),v). 
Jo Jm 1 Jo Jm 1 

Indeed, since u % is a solution, we have 

(3.32) I [ {Uu i -A{u i ){Vu\Vu i ),v) = -f [ ( J(u l )V(u l ), v), 
Jo Jm x Jo J Mi 

And the estimates (3.29), (3.30) holds true. So we have 

max II Vu l ||wfc-i,2 = max \\u\ + Vu l \\ W k-i,2 < C. 
te[o,T] te[o,T] 

when k > m + 1, by Sobolev, we know that for all t e [0, T] 

(3.33) W -> Vu in C°{M U N). 
and 

(3.34) Au* -> Am in L°°([0, T], L 2 (M 1? JV)). 
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The above convergence implies 

(3.35) lim / / {-Au i -A(u i )(Vu i ,Vu i ),v)=f [ (-Au - A(u)(Vu, Vu), v), 

l ->Wo J Mi JO J Mi 

and 

(3.36) lim f [ (J^Viu*)^) = [ f (J(u)V(u),v). 
^°°Jo Jm 1 Jo J Mi 

On the other hand, we have 

(3.37) lim/ T / (ul t ,v) = -[ T [ (u t ,v t )+[ ((u t (T),v(T)) — (u t (0),v(0))). 

l ^°°7o J Mi J0 J Mi J M x 

Now we can deduce from (3.32), (3.35), (3.36) and (3.37) that 
- ^ [ (u t ,v t )+ [ ((u t (T),v(T)) - (u t (0),v(0))) = 

JO J Mi J Mi 

[ [ (Au + A(u)(Vu,Vu),v) - I [ (J(u)V(u),v). 

Jo J Mi Jo J Mi 

This means u tt G L 2 ([0,T] x Mi,N), so we have proved (3.31), hence the theorem. 



□ 



Finally, we prove the uniqueness of the local solution. If u, v are two solutions to 
Cauchy problem (1.9), (1.10), we need to show u = v. Generally, one may consider the 
difference u — v between u and v. But in order to do the substraction, one needs to 
consider the embedding N ^ M. K . The following computation also relies on such an 
embedding. 

Theorem 3.10. Suppose u e W k ' 2 (M 1 , N),u x e W k ~ l ^(M u T U0 N) , where k > m for 
m > 2 and k = 2 for m — 1. Then the local solution to (1.9), (1.10) is unique in class 
W k < 2 . 

Proof. Assume u,v are two local solutions to (1.9), (1.10) satisfying 

u,v G L oo ([0,T],W k ' 2 (M 1 ,N)); u t ,v t G L°°([0, T],W k ~ 1,2 (M 1: TN)). 
Since we embed iV into a Euclidean space M. K , we can compute 

= [ (D t (u - v), D 2 (u - v)) - (D t (u - v), A(u - v)) 

J Mi 



Mi 

= [ (D t (u-v),(A(u)(Du,Du) - A(v)(Dv,Dv)) - (J(u)V(u) - J(v)V(v))) 

J Mi 

= / (u t , (A(u)(Dv,Dv) - A(v)(Dv,Dv))) - (v t , A(u)(Du, Du) - A(v)(Du, Du)} 

J Mi 

+ (ut - v t , -(J(u)V(u) - J(v)V(v))) 

< [ \A{u) - A(v)\((u t , \Dv\ 2 ) - (v t , \Du\ 2 )) +C f \u t - v t \\u - v\ 

JMi J Mi 

<C \u - v\\Du - Dv\(\Du\ 2 + \Dv\ 2 ) + C \u t - v t \\u - v\ 

J Mi J Mi 

< C\\Du-Dv\\ L 2 ■ (\\\u - v\(\Du\ 2 + \Dv\ 2 )\\ L 2 + ||«-u|| L2 ). 
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Hence we get 

(3.38) j t \\D(u - v)\\ L 2 < C(\\\u - v\(\Du\ 2 + \Dv\ 2 )\\ L , + || u - V || L2 ). 

If m < 3, we have k > 2. By Sobolev embedding W 2 ' 2 <->■ W 1 ' 6 , we get 

- V) \\ L 2 < C||« - v\\ L e(\\Du\\ 2 L e + ll^lUs) + C\\u - i;^ 

< C||Z)u - Dv\\v{\\Du\\wifl + ||Du|| w i, 2 ). 
If m > 3, we have Sobolev embedding W^ +1,2 W 1,2m . Thus 

(3 40) " W) " L2 " " (l|jDM|l ' 2m + H^ll^) + ^11" " V ^ L ' 2 

< C\\Du - Dv\\v(\\Du\\ wmfl + \\Dv\\ w w,2 + 1). 

From (3.39), (3.40) and Lemma 3.6, it follows that, if (uo,Ui) G W fc > 2 (Mi,JV) x 
W fc-1 ' 2 (M 1 , TN), there holds 

< C\\Du-bv\\ L 2. 

By Gronwall's inequality, we finally get 

\\b(u(t)-v(t))\\ L 2 < C\\b(u(0)-v(0))\\ L 2 = 0. 
Thus we complete the proof. □ 

Remark 3.1. We can also compute the difference between u and v intrinsically by using 
parallel translation. Mcgahahan [24] used this method to prove the continuous dependence 
of solutions to Schrodinger flow on initial data. Same method can by applied to prove 
continuous dependence of initial data to Cauchy problem (1.9), (1.10). 

Remark 3.2. We can also consider Schrodinger flow with potential, i.e. 

Ou 

(3.41) — = J(u)r(u) + J(u)VF(u), 

where F is a smooth function. Actually, we can prove the local existence of (3. 41) by the 
same method. 

4. Global existence in 1 + 1 dimension 

In this section, we follow the method in [31] to prove Theoreml.2. Note that when 
m — 1, m = 1 and k > 2 in Theorem 1.1. 

Proof of Theorem 1.2. According to Theorem 1.1, we already have a unique local solution 
u e L oo ([0,T),W 2 ' 2 (S 1 ,N)). Moreover, u satisfies the estimate (3.20). Now we need to 
derive a global estimate. Since u satisfies equation (1.9), i.e. 

(4.1) Uu = A(u)(Du, Du) - J{u)V{u). 

Applying a first order spatial derivative V to this equation, we get 

□(Vu) = V(A(u)(Du, Du)) - V(J(u)V(u)) 

= VA{u){Du, Du, Vu) + 2A{u){Vbu, Du) - J(u)WV(u) ■ Vu. 

But for the second fundamental form A, we have 

( Ut ,A(u)(;-)) = 0. 
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2dt 



\VDu\\ 2 L i = [ (□(V«),Vu t ) 

J Mi 



Thus 

(4.2) (Vu t , A(u){VDu, Du)) = (u t , VA(u)(VDu, Du, Vu)>. 
The above equality implies 

1 d 

||VDm||^ 2 = 

'Mi 

(VA(u){Du, Du, Vu) + 2A(w)(V£>w, Dm), V« t ) 

Mi 

-(j(M)vy(«)'Vti,Vw t ) 

= / (VA(u)(Z>u, Du, Vu), V« t ) - (J(u)VV(u) ■ Vu, Vu t ) 

J Mi 

+2 / (u t ,VA(u)(VDu,Du,Vu) 

J Mi 

(4.3) < C I \Du\ 3 \VDu\ + \Du\\VDu\. 

J Mi 

From Holder's inequality, 

(4.4) I \Du\ 3 \VDu\ < C\\Du\\ 3 L6 \\VDu\\ L 2. 

J Mi 

When m — 1, it follows from the classic Gagliardo-Nirenberg interpolation inequality and 
Kato's inequality that 

(4.5) \\Du\\ L 6 < \\VDu\\ a L2 \\Du\\ { lr a \ 
where 

- =a( --i )+( i- )-. 

i.e. a — \- Hence we arrive at a Gronwall-type inequality from (4.3), (4.4) and (4.5) 

^-\\VDu\\ 2 L 2 < \\V Du\\\4Du\\ 2 L 2 
at 

Combining this together with the energy inequality ||Du||| 2 < C, we obtain 

(4.6) \\VDu\\ 2 L 2 <C(t),Vt eR. 

Now we can derive the global existence from Theorem 1.1 and (4.6). Indeed, if this is 
not the case, assume the maximal existence time interval of u is [0, T). It follows from 
Lemma 3.6 that T only depends on the initial data, i.e. 



7-1,2 . 



r = T(||Du(0)||fli 

We may choose a small positive number e > 0, and consider the Cauchy problem (3.4) 
with initial data u(T — e). Then Theorem 1.1 guarantees the existence of another local 
solution u' G L°°([0,T'), W 2 < 2 (S\N)), where 

T' = T{\\bu{T-e)\\ H ^). 

Moreover, by the uniqueness Theorem 3.10, u and u' coincides on the overlapped time 
interval. Now, if we patch u, u' together, we get a solution to (1.9), (1.10) on the time 
interval [0, T — e + T"). The estimate (4.6) tells us that ||Z>u(£)||#i,2 is uniformly bounded 
for all t G [0,T). Consequently, if e is small enough, we have T — e + T' > T. This 
contradicts to the maximality of T. Hence, we must have T = oo. □ 
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Appendix A. Local existence of the approximation 

In this appendix, we use a fixed point argument to prove the local existence of the 
Cauchy problem of equation (3.5): 

(u tt - eAu t = F(u, u t ) 

' \«(0) = uo,«t(0) =ui 

where 

F(u, u t ) = Au - J(u)V(u) + A(u)(Vu + u t , Vw + u t ) - eT(u)(Au t ) 

and 

M e C7°°(Mi,JV),ui G C°°(M 1 ,TN) 
satisfy the following condition: 

u 1 (x)eT Uoix) N,\/xeM l . 

Consider the Banach spaces 

X = {v = ( Vl ,v 2 ) G C 3 (M U N) x C 2 (M u TN);v 2 (x) G T vl{x) N,Vx G Mi} 
with the norm 

IklU = lkl||c 3 (Mi) + ||^ 2 ||c* 2 (Mi) 

and 

y = c^M^JV) 

with the norm 

||/||y = ll/||cw 
We recall the expression of T(u)(Au t ) given by (3.6), i.e. 

(A.2) T(u)(Au t ) = A{u){Vu t , Vu) + dw{A{u){u t , V«)). 

From this equality, one can see that if (u,u t ) G X = C 3 x C 2 , then F(u,u t ) G C 1 . 
Therefore, F is a mapping from X into Y. In fact, we have 

Lemma A.l. F is a locally Lipschitz map from X to Y . 
Proof. For any v — (i>i, v 2 ), w — (wi,w 2 ) G X, we have 

\\F(v) - F(w)\\ Y < \\Avi - A Wl + J( Vl )V( Vl ) - JMVMWy 

+ \\A(v 1 )(Vv 1 + v 2 , Vui + v 2 ) - A(wi)(Vwi + u> 2 , Vwi + w 2 )||y 
+ e||T(ui)(Au2) -T(«;i)(A«; 2 )||y 
<I + II + III. 

Obviously, we have 

I < \\v — w\\x- 

For the second fundamental form, 

II < \\\A(v 1 )-A(w 1 )\\Vv 1 + v 2 \ 2 \\ Y 

+ ||A(iui)(|Vui - Vwi| + |v 2 - w 2 |)(|Vv + v 2 | + \Vw + w 2 \)\\ Y 
< C(\\v\\ 2 x + \\w\\ x )\\v - w\\ x . 
As for the third term, by a similar computation, we have 

/// < ep(ui)(Vu 2 , Vui) -i4(iui)(V«; 2 ,V«;i)||y 
+ e||div(^(ui)(u 2 , Vui) - A(iui)(«; 2 , Vwi))||y 
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Thus we obtain 

\\F(v) - F(w)\\ Y < C(l + \\v\\ 2 x + ||Hlx)lk - w\\ x , 
which means F is locally Lipschitz. □ 

It's well-known that there exists a heat kernel on compact manifold Ml, which we 
denote by H(x,y,t). We first fix u G X. Using the heat kernel, one can solve the linear 
parabolic equation 

f vt - eAv = F(u) 
(A ' 3) {,(0) = .i 

by 

v(x,t) = V(u)= H(x,y,et)u 1 (y)dy + / H(x,y,e(t - s))F(u{y))dydt. 
JMx Jo J Mi 

Then one can go on to solve an ordinary equation 

w t = *(«) 

w(0) = u . 



(A.4) 

The solution is given by 

w(t) = = / V(u)(s)ds + u . 
Jo 

Now we are ready to derive a fixed point argument. Fix 5 > 0, and set 
Z = {ue C([0,T],C 3 (Mi)) n C 1 ([0,T],C7 2 (M 1 )); 

(u,«t)|t=0 = («0,«l), ||(u(*),«i(t)) ~ ( M 0,Mi)|| X < 5} 

with the norm 

IMU = su p II («(*)>«*(*)) ll-x-- 

te[o,T] 

Lemma A. 2. $ : Z — > Z is a contraction ifT is sufficiently small. 

Proof. For any u,v E Z, we use the estimates of the heat kernel and Lemma A.l to get 

\mu) - $(v)\\ z < sup \\mu)-<S>(v),V(u)-V(v))\\ x 

te[o,T] 

< sup (|| r*(«(s))-*(u(s))tfe|| C 3 + ||*(«(*)) -*(u(*)||c») 



*G[0,T] JO 

3U 

= [0 



< sup / Ct~ a \\F{u,u t ) - F{v,v t )\\ Y ds 

*G[0,T] io 

< CT 1 "^ sup ||u(i) -u(t)||x, 



*G[0,T] 

where a G (0, 1) is a constant. Clearly if T is small, $ is a contraction of Z. □ 

Then by the Banach fixed point theorem, $ has a unique fixed point u £ Z, which is a 
local solution to equation (A.l). The regularity can be easily deduced from the property 
of the heat kernel. 
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